Integrable systems of quartic oscillators. 
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Abstract 

Several completely integrable, indeed solvable, Hamiltonian many-body prob- 
lems are exhibited, characterized by Newtonian equations of motion ("accelera- 
tion equal force"), with linear and cubic forces, in TV-dimensional space (N being 
an arbitrary positive integer, with special attention to N = 2, namely motions 
in a plane, and N = 3, namely motions in ordinary three-dimensional space). 
All the equations of motion are written in covariant form ( "N- vector equal N- 
vector"), entailing their rotational invariance. The corresponding Hamiltonians 
are of normal type, with the kinetic energy quadratic in the canonical momenta, 
and the potential energy quadratic and quartic in the canonical coordinates. 
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1 Introduction 



In a recent paper we introduced several integrable, indeed solvable, nonlinear 
oscillators. The basic idea was to start from the integrable, indeed solvable, 



and to transform it into a system of covariant, hence rotation-invariant, Newto- 
nian equations (" acceleration equal force" , with linear and cubic forces) for the 
time evolution of TV-vectors, by parameterizing appropriately the matrix U_ in 
terms of these TV-vectors. In the present paper we follow the same approach 1 
Q] HI, but we exploit certain additional, convenient parameterizations of matri- 
ces in terms of vectors, which we presented recently and which also transform 
the matrix evolution equation into covariant, hence rotation-invariant, New- 
tonian equations of motions for TV- vectors in TV-dimensional space. We thereby 
obtain - and exhibit in the following Section 2 - the Newtonian equations of 
motion of several new integrable, indeed solvable, systems of oscillators sub- 
jected to linear and cubic forces. As shown in Section 3, all these systems are 
Hamiltonian, with standard Hamiltonians of normal form: the sum of a stan- 
dard kinetic energy quadratic in the momenta and a potential energy depending 
(quadratically and quartically) only on the canonical coordinates. In view of 
the phenomenological relevance of such nonlinear oscillators, it is remarkable 
that several instances of such systems exist which are integrable indeed solvable. 
Investigation of their detailed dynamics is postponed to future papers. 

Notation: throughout this paper matrices are denoted by underlined char- 
acters, and vectors by superimposed arrows; their dimensionality will in each 
case be clear from the context. Superimposed dots indicate differentiations with 
respect to the independent variable ("time"), while dots sandwiched among two 
vectors denote the standard scalar product. As a rule, we write the equations of 
motion in the neatest form, forsaking any attempt to introduce additional con- 
stants by rescaling the independent variable and/or by redefining the dependent 
variables via linear transformations. 



2 Newtonian equations of motion 

A convenient way to obtain covariant ("vector equal vector") hence rotation- 
invariant Newtonian equations of motion is to insert in the matrix evolution 
equation Q the parameterizations introduced in jH] (hereafter the equation (n) 
of this paper |S] is identified here as (I.n)). Since the parameterization PI (1.3) 
was already exploited in pQ, we exhibit here only the new dynamical systems 
yielded by the new parameterization P2 (1.4). Let us emphasize that, since the 
matrix evolution equation Q is integrable, indeed solvable, all the dynamical 
systems obtained in this manner are as well integrable, indeed solvable. 
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Let us begin by writing down our basic result, obtained by inserting the 
formulas (1.4) into Q and by then performing some appropriate, merely nota- 
tional, changes: 
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Here, and throughout, the index n (as well of course as v) ranges from 1 to 
TV, and the index £ (as well of course as A) ranges from 1 to L, with TV and 
L two arbitrary positive integers; the 2 TV L time-dependent quantities ~r^ nt > = 
"^"^(t), f ( nt> = f ( nl \t) are TV- vectors (identifying the positions of moving 
point particles in TV-dimensional space); the dots sandwiched among TV— vectors 
denote the standard scalar product; and the (TV L) 2 quantities a^™ are arbitrary 
scalar constants. Because of the way they have been obtained, these 2NL 
Newtonian equations of motion, featuring linear and cubic forces, are integrable 
indeed solvable; and they are clearly covariant (TV- vector equal N- vector), hence 
rotation-invariant in TV-dimensional space. 

It would seem that this system always involves a number of TV- vectors which 
is an even multiple of the number of dimensions TV of the vector space. But this 
is not the case, due to the possibility of reductions, which may be enforced (with 
an appropriate choice of initial conditions) provided the constants oi? satisfy 
appropriate restrictions. For instance if these constants satisfy the symmetry 
restriction 



"■ex ~ "xe > 

it is clearly possible to impose the corresponding reduction 
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which obviously halves (from 2N L to N L) the number of evolving vectors, as 
well as the number of TV-vector Newtonian equations of motion, which clearly 
read then (see ©) 



(ni) 



N 



EE{" 

r/=l A=l 



>("A) 



(vn) 
l Xl 



>(„A) 



>{nty 



)]}■ 



(3c) 



Other (or additional) reductions of the Newtonian equations of motion J2J are 
moreover possible, in which some of the TV-vectors vanish identically, since 
clearly, provided for some value of the indices n and t the constants vanish 
unless v = n and A = I (a, » 
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5 n v Sex a>£ for v = 1, TV; A = 1, L with 
n, i fixed), then it is consistent, see (|2a.p. to set T^™^ = 0; and likewise, provided 



for some value of the indices n and £ the constants a 
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A = 1, L with n, £ fixed), 



and A = £ (<4a = &nv&t\<*t f° r v = !>•••> -^i 
then it is consistent, see (|2b)l . to set f ^ = 0. 

Let us also point out that, if all the constants vanish, so that the Newto- 
nian equations of motion only involve cubic forces, by applying an argument 
analogous to that used in [7] one concludes that all the nonsingular solutions of 
the following (to-deformed, complex) Newtonian equations of motion, 
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are completely periodic with period 2 7r / lu, provided m is a nonvanishing real 
(without loss of generality, positive) constant. Note that these equations of 
motion are complex (see the second term in their left-hand sides); but real 
equations can be obtained from these by introducing the real and imaginary 
parts of the N- vectors T > ( ne '> and f ( ne \ and of the "coupling constant" c (which 
has been introduced in @ via a trivial rescaling of the dependent variables), and 
by then considering separately the real and imaginary parts of these equations. 
There thus results a doubling (from 2 N L to 4 TV L) of the number of evolving 
real A^-vectors; but this number can be cut down via the same kind of reductions 
discussed above. No additional elaboration is reported below on these systems 
of nonlinear harmonic oscillators [?]• 

In the rest of this Section 2 we display explicitly, in the physically more 
interesting cases with N = 2 (evolutions in the plane) and N = 3 (evolutions in 
ordinary, three-dimensional, space), some simple examples involving few moving 
points. For the sake of completeness, and for the reader's enlightenment, we 
also display, in self-evident notation [S], the parameterizations that relate the 
equations of motion we report to the basic matrix evolution equation Q . 



2.1 Oscillators in the plane 

Let us consider the following parameterization of a 4 ® 4 matrix U (see Q) in 
terms of four 2— vectors (see (L6)): 
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as well as the following parameterization of the constant 4 (g> 4 matrix A: 



(a 











\ 




















a 







\o 


a( 12 ) 





a (22) _ 


a) 



(5d) 



Note that we are considering a case with N = 2, L = 1. 

We then obtain the following integrable, indeed solvable, system of four cou- 
pled linear plus cubic oscillators in the plane: 
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Further reductions to systems of less than four 2-vectors are easy to obtain by 
appropriate reductions. For instance by setting to zero the coefficient a' 21 ) and 
the 2-vector f ^ (and setting for notational convenience a' 11 ) = a, a' 12 ) = 
b, a' 22 ' = c and f W = "F^ 3 )) we get the following system of three coupled 
oscillators in the plane: 
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By setting moreover to zero "r^ 2 ) an d 6, and setting "F^ 1 ) = "r*, T^ 3 ) = f , we 
obtain the following system of two coupled oscillators in the plane: 
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And clearly the choice f — r reduces this system to the (trivially solvable) 
single equation 

7 > + a7 > = —^r r 2 . (9) 

A more general system of two coupled oscillators than (JHJ) obtains from JJjJ 
by setting a ( n ) = a, = a< 21 > = b, = c, ~?V> = 7«, 7^ = r>< 2 >: 
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And of course from this system the single equation © can be easily obtained 
via an additional, obvious, reduction. 



2.2 Oscillators in ordinary (3— dimensional) space 

Let us consider the following parameterization of a 6 (g> 6 matrix U (see (Q) in 
terms of six 3— vectors (see (1.7)): 
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as well as the following parameterization of the constant 6 ® 6 matrix ^4: 
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Note that we are considering a case with N — 3, L = 1. 

We then obtain the following integrable, indeed solvable, system of six cou- 
pled linear plus cubic oscillators in ordinary (3— dimensional) space: 
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It is plain that if a^ n ^ — a^ nv ^ one can set f (") = "r* 1 ™), obtaining thereby 
the following system of i/iree coupled linear plus cubic oscillators: 
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Reductions of the system (|12|l to a number of oscillators smaller than six can 
be simply obtained by setting to zero some vectors and, for consistency, some 
of the constants a (rll/ ). For instance, by setting a^ 31 ' = a (13 ) = a^ 32 ) = a' 23 ) = 
(and, for notational convenience, a^ 11 ' = a, a' 21 - 1 = 6, a^ 12 ) = c, a^ 22 ) = d) 
and T*( 3 ) = f ^ — 0, one gets the following system of four coupled linear plus 
cubic oscillators in ordinary (3-dimensional) space: 
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Note that, denoting by tt the plane identified by the requirement that the 
two vectors "r^ 1 ) and "r^ 2 ) lie initially (namely, at t = 0) in it, if the initial 

_^(l) • (2) 

velocities 7 (0) and "r* (0) also lie in the same plane 7r, then the two vectors 
7^ 1 -'(t) and 7( 2 )(t) remain in this same plane 7r throughout their time evolution 
- independently of the behavior of the two vectors f W and f ^ , which of 
course might as well also always remain on some (possibly different) plane, say 
7r, if their initial velocities lie in the same plane defined by their initial positions. 

Further reductions to two, or just one, oscillators yield equations (for 3- 
vectors) analogous to those written above for 2- vectors (see (JHJ), © and (|TU)l % ). 
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3 Hamiltonians 



The matrix evolution equation Q is Hamiltonian: it obtains, in the standard 
manner, from the Hamiltonian 



H{P,U_) = trace 
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where the canonical coordinates are the components of the matrix U_ and the 
corresponding canonical momenta are the components of the matrix P. (Indeed, 
the integrability of this matrix evolution equation, ijTjl. is due to its being a 
reduction HU |2| @] of the Non-Abelian Toda Lattice [S] the Hamiltonian 
structure of which was already exploited in This clearly entails that all the 
Newtonian equations of motion reported in this paper (except are as well 
Hamiltonian. 

We display here only the Hamiltonian function that, via the standard Hamil- 
tonian equations, yields the basic Newtonian equations of motion It reads: 
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In this Hamiltonian, the canonical coordinates are the 2 N 2 L components of 
the 2NL iV-vectors r'™'',^" 1 ), n = 1,2,..,N, I = 1,2,..,L, and the corre- 
sponding canonical momenta are the 2 iV 2 L components of the 2 N L N- vectors 
pHn^ p n = l,2, .., N, 1=1, 2, .., L. Note that this Hamiltonian entails 
that 
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Also note that this Hamiltonian is of normal type, namely it features a kinetic 
energy term depending quadratically only on the momenta and a potential en- 
ergy term depending (quadratically and quartically) only on the coordinates. 



4 Final remarks 

As we already emphasized above, the Newtonian equations of motion considered 
in this paper are all covariant hence rotation-invariant, this being a rather 
essential condition to interpret them as describing a many-body problem. They 
are, however, not invariant under translations; but variants of them that do 
possess this additional property, without forsaking the property to be integrable 
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indeed solvable, can be easily manufactured (for a technique to do so see @]). 
This generalization, as well as detailed explorations of the actual behavior of 
the solutions of these models, are left as tasks for the future; as well as the 
application of these findings to phenomenologically interesting situations (in 
physics or elsewhere). 

A more challenging task for the future - hopefully appealing to physicists 
and to mathematicians - will be the study of the quantal counterparts of the 
models discussed in this paper. 
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